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Abstract - Using a sum rule approach we investigate the dipole oscillation of a spin-orbit cou- 
pled Bose-Einstein condensate confined in a harmonic trap. The crucial role played by the spin 
polarizability of the gas is pointed out. We show that the lowest dipole frequency exhibits a char- 
acteristic jump at the transition between the stripe and spin-polarized phase. Near the second 
order transition between the spin-polarized and the single minimum phase the lowest frequency is 
vanishingly small for large condensates, reflecting the divergent behavior of the spin polarizability. 
We compare our results with recent experimental measurements as well as with the predictions of 
effective mass approximation. 



Synthetic gauge fields represent a rapidly developing 
direction of research in ultra cold atomic physics both 
from the experimental jlTHO] and theoretical perspective 
[11H26] . They give rise to the occurrence of new quantum 
phases exhibiting unique magnetic features. In the pres- 
ence of spin-orbit coupling also the collective oscillations 
exhibit challenging features which have been already the 
object of theoretical studies [27H31] . as well as first exper- 
imental measurements [5] . In particular the experiment of 
[5] has shown that the center-of-mass oscillation of a har- 
monically trapped Bose-Einstein condensate can be deeply 
affected by the coupling with the spin degree of freedom. 
In the same paper the problem was investigated theoreti- 
cally using a variational approach within time dependent 
Gross-Pitaevskii theory. 

The purpose of the present work is to study the be- 
havior of the center-of-mass oscillation employing a sum 
rule approach. This method is known to point out the 
role of conservation laws and to reduce the calculation of 
the dynamical properties of the many-body system to the 
knowledge of a few key parameters relative to the ground 
state. In the case of the center-of-mass oscillation we will 
show that, due to the spin-orbit coupling, a key role is 
played by the spin polarizability of the gas which is re- 
sponsible for important deviations of the dipole frequency 
from the harmonic oscillator value. 

In this work we consider a spin 1/2 interacting gas 
characterized by the Hamiltonian (for simplicity we set 
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where i = 1 , ■ • ■ , N is the particle index, while a and /3 
are the spin indices (t, 4- = ±) characterizing the two spin 
states. The single-particle Hamiltonian ho is defined by 



/in = r 



{Px-k (J z ) +p 2 ± 



o s 

2°x + ^z + V cxt (r) (2) 



and is characterized by equal contributions of Rashba |32j 
and Dresselhaus [35] spin-orbit couplings and a uniform 
magnetic field in the (x, z)-plane, with the Raman cou- 
pling constant accounting for the transition between the 
two spin states, fco the momentum transfer of the two Ra- 
man lasers, 5 the energy difference between the two single- 
particle spin states, and V ex t = (oj 2 x 2 +uj 2 y 2 +uj 2 z 2 )/2 the 
external trapping potential, thereby chosen of harmonic 
type. The spin up and down density operators entering 
Eq. Q are defined by n±(r) = (1/2) £. (1 ± a z ,i) S(r-r t ) 
while g a fj = ATta a p are the relevant coupling constants in 
the different spin channels, with a a p the corresponding s- 
wavc scattering lengths. Finally cr^, with k = x,y,z, are 
the usual 2x2 Pauli matrices. The Hamiltonian ([T]), ([2} 
has been already implemented experimentally [US] . It has 
been recently employed to describe a variety of non trivial 
quantum phases in Bose-Einstein condensates [T8|[T9j . 
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In the absence of spin-orbit coupling (fco = 0) the center- 
of-mass oscillation of the trapped gas along the ^-direction 
is exactly excited by the operator X = J2i x i an d corre- 
sponds to the usual sloshing mode with frequency uj x . In 
order to investigate the effect of spin-orbit coupling it is 
useful to introduce the following moments of the excitation 
strengths of the operator F at zero temperature 0: 



{F)=Y,{E n -E ) k \(0\F\n){ 



(3) 



where |0) and |n) arc, respectively the ground state 
and the excited states of H with E n — Eq the corre- 
sponding excitation energies, while | (0 1 ^ | ^) | 2 is the F- 
strength relative to the state n). Some relevant sum 
rules can be easily calculated employing the closure re- 
lation and the commutation rules involving the Hamil- 
tonian of the system. For example the most famous 
energy-weighted moment for the dipole operator F = X 
takes the model independent value (also called /-sum rule) 
mi{X) = (1/2) (0| [X, [H,X]] |0) = N/2 where N is the 
total number of atoms. Notice that this sum rule is not 
affected by the spin terms in the Hamiltonian, despite the 
fact that the commutator of H with X explicitly depends 
on the spin-orbit coupling: 



[H,X) = -i(P x -k X 2 



(4) 



Here P x = X^Pz.i i s the total momentum of the gas along 
x and E z = J2i a z,i is the total spin operator along z. 
Equation ((4]) actually reflects the fact that the equation of 
continuity (and hence in our case the dynamic behavior of 
the center-of-mass coordinate) is deeply influenced by the 
coupling with the spin variable. The above commutation 
rule will be later used to exploit the importance of such a 
coupling in the evaluation of the m_3(A) sum rule. 

Another important sum rule is the inverse energy- 
weighted sum rule, also called dipole polarizability. In the 
presence of harmonic trapping this sum rule can be calcu- 
lated in a straightforward way using the commutation re- 
lation [H, P x ] = iui x X and the closure relation. One then 
finds the result rri-i(X) = N/ {2lo x ), thereby showing that 
both the energy- weighted and the inverse energy- weighted 
sum rules relative to the dipole operator X are insensitive 
to the presence of the spin terms in the single-particle 
Hamiltonian as well as to the two-body interaction. 
This does not mean, however, that the dipole dynamic 
structure factor is not affected by the spin-orbit coupling. 
This effect is accounted for by another sum rule, partic- 
ularly sensitive to the low energy region of the excitation 
spectrum: the inverse cubic energy-weighted sum rule for 
which we find the exact result 
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where we have introduced the spin polarizability per par- 
ticle x = 2to_i(E 2 )/A relative to the z-th direction of 



x At finite temperature, the moments should include the 
proper Boltzmann factors 1341 



the spin operator. In order to derive result ([5]) we have 
used the recurrence relation m_ 3 (A) = m-i(P x )/u x fol- 
lowing from the commutation relation for [H, P x ], as well 
as the most relevant commutation rule (UJ) for [H, X] . It 
is worth mentioning that the above results for the sum 
rules mi(X), m_i(X) and m-%{X) hold exactly for the 
Hamiltonian (TTJ, including the interaction terms. Their 
validity is not restricted to the mean-field approximation 
and is ensured for both Bose and Fermi statistics, at zero 
as well as at finite temperature. In particular the sum rule 
ms(X), being sensitive to the magnetic susceptibility, is 
expected to exhibit a non trivial temperature dependence 
across the BEC transition. 

Equation ([5]) exploits the crucial role played by the spin- 
orbit coupling proportional to fco. The effect is particu- 
larly important when the spin polarizability takes a large 
value. A large increase of \ is associated with the occur- 
rence of a dipole soft mode as can be inferred by taking 
the ratio between the inverse and cubic inverse energy- 
weighted sum rules m_i(X) and m_ 3 (A), yielding the 
rigorous upper bound uj x / y/l + k^x to the lowest dipole 
excitation energy. The effect of the coupling between the 
center-of-mass oscillation and the spin degree of freedom 
is further revealed by making the ansatz F = P x + ?7fcoE z 
for the optimized operator exciting the dipole oscillation 
and minimizing the collective frequency fixed by the ratio 
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with respect to variations of the real parameter rj. In de- 
riving the last equality we have explicitly used the sum rule 
result mi(S,) = (1/2) (0| [E z , [H,E Z ]] |0> = -NSl(a x ) for 
the energy-weighted moment relative to the spin excita- 
tion operator E z . The choice rj = in Eq. ^ corresponds 
to the estimate m_i(X)/m_ 3 (A). In the opposite 77 ^> 1 
limit Eq. ([6]) instead coincides with mi(E 2 )/m_i(E z ). Re- 
sult ([6|), based on the linear response formalism of sum 
rules, provides a rigorous upper bound to the dipole fre- 
quency in the regime of small amplitude oscillations. 

Let us now discuss the behavior of the physical quanti- 
ties x an( i ( a x) entering the expression for the dipole fre- 
quency. Their actual behavior depends on the quantum 
phase characterizing the ground state of the many-body 
system. These phases were investigated in details in [19] 
for a spin 1/2 interacting Bosc-Einstcin condensate em- 
ploying the ansatz 
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(7) 

for the order parameter. The variational calculation, ap- 
plied to a uniform configuration, yields the conditions 
k± = ±fci and 8± = 6 = arccos(fci/fco)/2 and predicts 
the occurrence of three quantum phases, depending on 
the value of the relevant parameters fco, and g a p (here 
and in the following we consider a spin symmetric Hamil- 
tonian with <7-|-|- = <7_|4 = g and 5 = 0). These arc the stripe 
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phase (also called phase I) where Bose-Einstein condensa- 
tion takes place in a combination of plane waves with op- 
posite momenta ±fci, provided g > g^±, the spin-polarized 
phase II and the single minimum (ki = 0) phase III. A nat- 
ural generalization of the variational wave function can be 
also used to calculate the dynamic properties of uniform 
matter, by solving the Bogoliubov equations and deter- 
mining the relevant elementary excitations of the system, 
the momentum distribution and the quantum depletion of 
the condensate [35] . 

The static spin polarizability of the spinor Bose-Einstein 
condensate at zero temperature can be directly calculated 
by expanding the values of k± and 0± around the equilib- 
rium values and minimizing the energy of the system in the 
presence of an external static field proportional to a z . For 
densities smaller than the critical value — k 2 /(2"/g), 
with 7 = (g — g^i)/ (g + g^]), all the three phases discussed 
above are available [T5]. In the stripe phase, holding for 
small Raman frequency SI, one finds the result 



n 2 - 4fc4 



(d + 2G 2 )o 2 -8G 2 i4 



(8) 



where we have defined G\ = n(g + g^i) /4 and G2 = 
n{g — g^\) /4 and, for simplicity, we have considered 
the weak coupling limit characterized by the condi- 
tion G\, G 2 <C k 2 . The spin polarizability xi di- 
verges as one approaches the critical frequency fit 1 " 11 ) = 
2kQ^/2j] (1 + 27) providing the transition to the spin- 
polarized phase. However Eq. (JS| is valid only in the low- 
density limit and inclusion of higher-order terms makes the 
value of x finite at the transition. In the spin-polarized 
phase II the spin polarizability is given by 



o 2 



(k 2 -2G 2 ) A{k 2 -2G 2 y 



0- 



(9) 



It takes a larger and larger value as one approaches the 
transition point to the zero momentum phase III occur- 
ring at the frequency O^ 11-111 ) = 2{k 2 — 2G 2 ). In the zero 
momentum phase III the spin polarizability takes instead 
the value 

Xln(fi) = o oJ ( 10 ) 
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2(k 2 
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exhibiting a divergent behavior when one approaches the 
transition point from above and vanishing in the limit of 
large O. Results © and (JTUJ), whose validity is not limited 
to the weak coupling regime, explicitly reveal the second 
order nature of the phase transition, the values of x dif- 
fering by a factor 2 when one approaches the transition 
from above or below. They also show that the relevant 
interaction parameter in both phases is the spin density 
dependent parameter G 2 ■ The behavior of the spin polar- 
izability as a function of the Raman coupling O is shown 
in FigfTJ (a) for a choice of parameters which emphasizes 
the deviations with respect to the weak coupling result. In 
particular the transition point between the phases II and 




Fig. 1: (Color online) (a) Spin polarizability \ as a func- 
tion of Q, calculated in uniform matter (blue solid lines) and 
in the harmonic trap (red dashed lines). (b) The corre- 
sponding lowest mode frequency w with (a x ) and \ calcu- 
lated in uniform matter (blue solid lines) and in the har- 
monic trap (red dashed lines). The black dotted lines cor- 
respond to the prediction cu x /^/l + k^x ( see text). The pa- 
rameters: kg — 2-7T x 320 Hz, uj x = 2-k x 20 Hz, density in 
the center of the trap n ~ 2.6 x 10 13 cm -3 , scattering length 
a tt — a ll = 100 as, a^i — 60 as, where as is Bohr radius, 
corresponding to Gi/fco — 0.257 and Gijkk — 0.064. 



Ill is predicted to be located at a value ~ 15 % smaller 
than the value 2&q (the weak coupling result). In the 
same figure we report the value of the spin polarizability 
calculated within Gross-Pitaevskii theory in the presence 
of ID harmonic trapping [19) . The comparison with the 
uniform matter prediction, corresponding to the value of 
the density calculated in the center of the trap, is rather 
good reflecting the usefulness of the uniform matter cal- 
culation of x- 

The behavior of the average transverse spin polariza- 
tion (<t x ) was calculated in |19j . In the stripe and in 
the spin-polarized phase the uniform matter calculation 
yields, respectively, (a x ) = (0/2) / (fcg + Gi) and (a x ) = 
(0/2) / (ky — 2G2) , while in the single minimum phase 
one has (a x ) = — 1. 

We are now in a position to estimate the dipole fre- 
quency by minimizing Eq. ([6]) with respect to the param- 
eter 77. The variational procedure actually provides two 
solutions. The upper solution is physically meaningful 
only for very small values of O where it approaches the fre- 
quency u> x of the center-of-mass sloshing mode. For higher 
O the upper solution takes large values and the coupling 
with other modes, not accounted for by our ansatz for the 
excitation operator F, becomes important. The results 
for the lowest dipole solution are reported in FigQ] (b). 
They reveal the important deviations from the oscillation 
frequency ui x caused by the spin-orbit and Raman cou- 
plings for all values of O. In the same figure we also show 
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the prediction ui x / \/l + k^x for the dipole frequency, ob- 
tained by setting 17 = in {BJ. This turns out to be an 
excellent estimate except for very small values of the Ra- 
man coupling. Actually, in the 57 <C uj x limit the value 
of r\ minimizing Eq. ^ is no longer small and the mode 
turns out to be a pure spin oscillation, its frequency van- 
ishing linearly with 57. In this limit the mode does not 
exhibit any significant coupling with the center-of-mass 
oscillation. At the transition between the phases I and II 
the lowest dipole frequency exhibits a sudden jump and 
then starts decreasing for larger values of 57. In the ther- 
modynamic limit, it vanishes at the transition between the 
phases II and III as a consequence of the divergent behav- 
ior of x and, above the transition, it increases to reach 
asymptotically the oscillator value ui x at large 57. 

In FigJ5] we show the spin polarizability and the fre- 
quency of the lowest dipole mode calculated in the experi- 
mental conditions of [8] . They correspond to a very small 
value of 7 (~ 10~ 3 ), thereby causing the compression of 
the stripe phase into a narrow region at small values of 57. 
Furthermore these conditions correspond to a very small 
value of G2. As a consequence, in the wide range of Ra- 
man coupling 57 3> to x Eqs. (|6|), ([9]) and (p~0|) yield the 
useful results 

c4 = col (1 - ^ 2 /4fc 4 ) , w? n = u*(l-2*g/fi) (11) 

for the dipole frequency in the phases II and III. FigJ5] (a) 
shows in a clear way the divergent behavior of the spin po- 
larizability at the transition between the two phases which 
is responsible for the quenching of the dipole frequency. 
The GP simulations are practically indistinguishable from 
the calculations of x m the uniform matter (and hence of 
lv) based on uniform matter ingredients. The comparison 
with the experimental data for the dipole frequencies, as 
shown in Figj2](b), is good far from the transition point at 
57 = 2&Q, while near the transition nonlinear effects play a 
major role as discussed in [5] (see also discussion below). 

The combined spin-orbit nature of the lowest dipole 
mode is also nicely revealed by the relative amplitudes 
of the oscillating values of the center-of-mass position 
(Ax), momentum (Ap) and spin polarization (Ay). These 
amplitudes can be calculated in the present approach 
by writing the many body oscillating wave function as 
= e Mt)F e mo | ) where F = Px+ vkQ z z i s the 

excitation operator and D is defined by the commuta- 
tion relation [H, D] = F, while a and /3 are parame- 
ters whose time dependence can be obtained through a 
variational Lagrange procedure. The time dependence of 
the relevant quantities (X), (Px) and (£ z ) is then eas- 
ily calculated by expanding the wave function up to first 
order in a and f3. For 57 <C lj x , as discussed above, 
the lowest frequency mode is mainly a spin oscillation 
(large rf) and one finds that the center-of-mass position 
is basically at rest (Ax ~ 0), while Ap/kg ~ Ay,. For 
larger values of 57 the lowest frequency is instead well ap- 
proximated by the choice r\ = 0, and the relationships 
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Fig. 2: (Color online) (a) Spin polarizability x an d (b) lowest 
mode frequency cu as functions of 57 (see Fig[TJfor details) with 
different choice of parameters: k$ = 2tt x 4.42 kHz, ui x — 2ir x 
45 Hz, density in the center of the trap n ~ 1.37 x 10 14 cm -3 , 
scattering length aft = 044 = 101.20 as, 0^-4. = 100.99 as, 
where as is Bohr radius, and the atomic mass of 87 Rb. The 
black circles are the experimental data of [H] . The black arrow 
in (a) indicates the transition between the phase I and II. 

between the spin, center-of-mass and momentum ampli- 
tudes take the useful form Ay = AxkoU) x x/\/l + k^x 
and Ap/ko = Ay (l + k^x) I (&oXj- The connection be- 
tween the momentum and spin amplitudes has been al- 
ready pointed out in [8]. The above relationships show 
that, near the transition between the phases II and III, 
where the spin polarizability diverges, the amplitude Ay 
of the spin oscillation can easily take large values, thereby 
emphasizing the role of nonlinear effects. These effects are 
likely at the origin of the finite values of the dipole frequen- 
cies observed at the transition (see Fig. 2 and discussion in 
[5]). It is finally worth mentioning that, despite its strong 
spin nature, the lowest frequency mode exhausts almost 
completely the dipole polarizability sum rule m_i(X), ex- 
cept in the Q -C u x region. As a consequence it can be 
easily excited by displacing the trapping potential. 

In the last part of the work we show that a useful insight 
on the lowest dipole oscillation can be also obtained by 
investigating the dynamic behavior employing directly the 
lower branch 

<P) = I (Pl+Pl + fco) - + S! 2 (12) 

of the excitation spectrum of uniform matter. The new 
single particle Hamiltonian can be used to solve the classi- 
cal Hamilton- Jacobi equations in both the linear and non- 
linear regimes [8]. It can be also used to calculate the 
ratio mi(X)/m—i(X) between the energy-weighted and 
inverse energy-weighted sum rules relative to the dipole 
operator X. While the result N/(2ia x ) for the dipole po- 
larizability continues to hold, the energy-weighted moment 
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is deeply affected by the new Hamiltonian reflecting the 
projection procedure accounted for by Eq. (fT^j) . By car- 
rying out explicitly the double commutator one finds the 
result us 2 = u;^{<9^e(p)). By replacing the value of p x with [10 
the value where e(p) is minimum, i.e., where Bosc-Einstein 
condensation takes place in uniform matter (effective mass 
approximation) , one recovers the results predicted by 
the sum rule approach in the weak coupling limit in the 
phases II and III respectively. A better estimate of the fre- 
quency at the transition point f2 = 2/cq can be obtained by 
evaluating the integral J dp n(p) dp x e(p) with n(p) given 
by the momentum distribution of the Gross-Pitaevskii so- 
lution in the presence of the harmonic trap. With the 
choice of parameters of Fig[2]we find U) m - m ~ 0.03 u x . The 
use of lowest branch dispersion (TT^|) turns out to be instead 
unsuitable to investigate the lowest dipole frequency in the 
stripe phase. 

In conclusion we have developed a sum rule description 
of the dipole oscillation which explicitly reveals the cru- 
cial role played by the spin polarizability in the presence 
of spin-orbit coupling. The method can be also applied 
to Fermi gases as well as generalized to study the cou- 
pling between the dipole and other collective oscillations 
recently pointed out in [31] . 

Note added after submission: the divergent behavior of 
the magnetic susceptibility predicted in the present work 
has been recently observed experimentally in [36) through 
the analysis of the relative spin and momentum ampli- 
tudes. 
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